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Abstract: We consider the effect of a periodic perturbation with frequency to on the 
^ ' holographic = 4 plasma represented by the planar AdS black hole. The response of 

o 
in 

. numbers can be found by solving wave equations in the AdS black hole background with 

, infalling boundary conditions on the horizon in an analogous way as in the calculation 

\ of quasinormal modes. The complex momentum eigenvalues have an interpretation as 

\ poles of the retarded Green's functions, where the inverse of the imaginary part gives an 

' absorption length A. At zero frequency we obtain the screening length for a static field. 

These are directly related to the glueball masses in the dimensionally reduced theory. We 

^ ' also point out that the longest screening length corresponds to an operator with non- 

^ . vanishing R-charge and thus does not have an interpretation as a QCD3 glueball. 
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1. Introduction 

The AdS/CFT correspondence is a concrete realization of the the idea that the large- A'^ 
limit of non-Abelian gauge theories can be described by a dual string theory [1]. In the 
large 't Hooft coupling regime the dual theory admits a description in terms of gravitational 
fields over a weakly curved background. More precisely the AdS/CFT correspondence pro- 
poses an exact duality between M = 4 super Yang-Mills in four dimensions and type IIB 
superstrings in AdSs x S^ . This theory is conformal, but at finite temperature confor- 
mal symmetry is broken and the theory is in a deconfining (or plasma) phase. The dual 
description corresponds to a black hole geometry with flat horizon [2]. 

Heavy-ion collisions at RHIC [3] and lattice simulations [4] indicate that QCD actually 
stays strongly coupled above the deconfinement transition up to temperatures T ~ 2Tc. 
Therefore, it is of great interest to develop non-perturbative tools that can describe the 
properties of the strongly coupled plasma. Lattice simulations are good to describe ther- 
modynamical properties, but out-of-equilibrium processes are much harder to analyze. In 
this context, the AdS/CFT correspondence could provide a better framework to derive 
some plasma properties using analytic methods. Although M = 4 is far from QCD, some 
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qualitative properties of the plasma seem to be quite similar, and AdS/CFT computations 
of the shear viscosity [5], the energy loss rate of a heavy quark [6] or the jet quenching 
parameter using light-like Wilson loops [7] show good agreement with experimental data. 

Many of the properties of the plasma can be studied using linear response theory. In 
this approximation, small perturbations that do not change significantly the state of the 
plasma are introduced. The system then tries to restore thermal equilibrium. That involves 
dissipation if the perturbations are localized in time or absorption if they are localized in 
space. ^ We would like to address the latter in this work. The absorption is directly related 
to spatial correlations in the equilibrium state. At high temperatures the system is in a 
very disordered phase, so measures made in different parts of the plasma give uncorrelated 
results. For the same reason, a small perturbation cannot travel too far in the plasma 
before being washed out by thermal fluctuations. How far this can be depends on the 
details of the plasma, but in general we expect that the characteristic absorption lengths 
decrease as the temperature increases. 

In the gravity dual picture the absorptive properties of the plasma rely on the presence 
of a horizon. Small classical perturbations end up falling into the horizon, either after a 
finite time or after travelling a finite distance.'^ The first is described by complex values for 
the eigenfrequency, the quasinormal modes, and the second by complex momentum values. 
Both are intimately related, they correspond to solutions of the linearized equations of 
motion. They also satisfy the same boundary conditions, Dirichlet at the AdS boundary^ 
and infalling at the horizon. The difference is that quasinormal modes decay exponentially 
in time while complex momenta describe the decay along the direction of propagation. 
The choice of boundary conditions restricts the possible values of complex frequency or 
momentum to a discrete set. In the dual gauge theory we can interpret them as inverse 
relaxation times r or inverse absorption lengths A of the plasma. The relaxation time 
depends on the (real-valued) wave number k whereas the absorption length depends on the 
(real- valued) frequency uj, i.e. A = X{lo). In the gravity theory we therefore fix a frequency, 
impose infalling boundary conditions on the horizon and search then for a solution of the 
boundary condition at infinity in the complexified momentum plane. In this way we can 
compute the frequency dependence of the absorption lengths. 

Quasinormal modes have been much studied in the context of black holes in flat space- 
time (see [9] for a review) and in the AdS/CFT correspondence [8,10-15]. Complex mo- 
menta have been studied for horizons of compact spatial geometry, where they correspond 
to Regge poles [16] of the black hole S-matrix. In the AdS/CFT correspondence, the zero 
frequency limit of complex momenta gives the glueball masses of QCD3, as computed in 
refs. [17-19]. While this work was in progress the interpretation of glueball masses as 
correlation lengths was also emphasized in [20]. 

The content of the paper is the following. In section |2| we explain in detail the relation 
between complex momenta and absorption lengths. We show that they arise as the poles 

""^ There can also be diffusion effects if conserved charges are involved. 
^In AdS space the curvature acts effectively as a box, so they cannot escape to infinity. 
^Actually, the condition is that they should be normalizable modes, so in the gauge theory they corre- 
spond to states and not to the insertion of sources or couplings (c.f. [8]). 
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of the retarded Green's function and give an argument based on stability considerations 
showing that the poles have to lie in the first and third quadrants of the complex momentum 
plane. In section ^, we compute the frequency dependence of the largest correlation lengths 
for scalar operators of conformal dimension A = 4, global currents, and the transverse and 
shear channels of the stress-energy tensor, respectively. We also show the relation with 
QCD3 glueball masses. The paper ends with a summary of our results and some outlook 
to future possible investigations in section In appendix ^ we comment on the form of 
the effective potentials that arise in rewriting the wave equations on AdS in the form of a 
Schrodinger equation and in appendix ^ we show how to avoid the "false frequencies" that 
arise in the numerical algorithm based on the Heun equation [13]. 

2. Absorption lengths in AdS/CFT 

In ref. [22] the authors gave a prescription to compute retarded two-point Green's functions 
in the context of the AdS/CFT correspondence with Lorentzian signature. It was empha- 
sized that retarded propagators correspond to imposing an infalling boundary condition 
at the horizon for the fields on the gravity side. On the other hand, infalling boundary 
conditions are also the constitutive ingredient for the calculation of the quasinormal fre- 
quencies of black holes in anti de Sitter space [10].^ The authors of [11] observed that 
the quasinormal frequencies of BTZ black holes coincide with the poles of the retarded 
two-point functions in the dual two dimensional conformal field theory. In [8, 13] it was 
shown that this observation extends generally to the Lorentzian AdS/CFT correspondence, 
i.e. quasinormal frequencies in AdS can be interpreted as the poles of retarded Green's 
functions in the dual field theory. 

Let us remember the interpretation of the poles of the retarded Green's function Gr. 
The response in the field (p of the system under consideration is obtained by the convolution 
of the perturbation represented by the source j{t,x) with the retarded Green's function 

^{t, ^)-=- j dr d^ Gn{t - r, X - ^) i(T, ^) . (2.1) 

If one chooses a perturbation localized in time, figuratively speaking one "hits" the plasma 
once at time i = 0, the perturbation is given by j(t,x) = (5(t) exp(iqx).^ Considering 
the Fourier transform of the retarded propagator and performing r and ^ integrations one 
arrives at 

0(t,x) = j GR(i/,q)e— * . (2.2) 

One can now make the analytical continuation to the complex z/-plane and use Cauchy's 
theorem. For t > we form a closed contour with a semicircle at infinity on the lower-half 
z^-plane, whereas for t < we would close it in the upper-half i/-plane. One obtains 

(/)(t,x) =isign(t)e*'i^ V e"**^"* ResGR(z/,q) , (2.3) 



*For a more general review of quasinormal modes see [9]. 

^The a;-dependence is that of a plane wave; however a general dependence can be constructed by super- 



positions of plane waves. 
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At this point we assume that the retarded Green's function is analytic in the upper half 
of the frequency plane and that its only singularities are single poles in the lower half 
plane. This is indeed the analytic structure that appears in the Lorentzian AdS/CFT 
correspondence at finite temperature [8]. In general, the analytic structure of retarded 
two-point functions is of course more complicated and involves also branch cuts. The 
authors of [21] computed the retarded two-point function of tr (-F^) at weak coupling and 
found a tower of branch cuts with branch points located at u; it g = —lAimT. In this paper 
we will only consider the strict large and strong coupling limit. Therefore, the response of 
the system to a perturbation localized in time is determined by the sum over the residues of 
Gr at the poles. In the holographic gauge theory these poles are precisely the quasinormal 
frequencies of the perturbation on AdS space subject to the infalling boundary condition. 

Instabilities, i.e. exponentially growing modes, appear as quasinormal frequencies with 
positive imaginary part. This is consistent with the interpretation as retarded Green's 
function, where singularities in the upper half plane would correspond to tachyonic modes 
travelling backwards in time. A typical arrangement of quasinormal frequencies as they 
appear in the analysis of small perturbations of asymptotically AdS black hole spacetimes 
is depicted in figure ||. 



Im u 



Re u 



Figure 1: The relevant integration contour for the poles in the z^-plane. All the poles are in the 
lower-half plane, corresponding to the interpretation of quasinormal modes as the poles of a retarded 
Green's function. 



Let us choose now another kind of perturbation. This time we will pick a periodic 
perturbation localized in space, i.e. we switch the roles of time and one space coordinate 
and assume a source of the form j{t, x) = 5{x) exp[—i{ujt — k_|_x^)]. We compute the effect 
of such a perturbation again in linear response theory. Doing the Fourier transform of the 
retarded propagator and performing and r integrations one finds 

^^^'""^ = -(i) / ^nioJ^k^,q) e^"" , (2.4) 

This is the response of the system to a periodic perturbation with frequency uj that is 
localized in the j;-direction and has the form of a plane wave in the perpendicular directions 
x_L. We have assumed that the perturbation has started far in the past such that all 
transient oscillations have already vanished and the system has reached a stationary state. 
In the following we will also assume that the perturbation is not further modulated in the 
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x_|_-directions, i.e. we set k_|_ = 0. Now one can use again Cauchy's theorem, closing in the 
upper or lower-half planes for x > and x < 0, respectively. The result is 



(p{t, x) = —i sign(x) e' 



-iujt 



E 

(j„;poles 



e^"?"^ Res Gr(w, 



q=qn 



(2.5) 



Again we have assumed that the only singularities of Gr(u;, q) are poles in the complexified 
momentum plane. By symmetry considerations {x —x) it is clear that q = qn = 
Qn + is a pole then also q = —qn has to be a pole. We would like the poles in the 
upper-half to lie in the first quadrant and those in the lower-half in the third quadrant. 
With such and arrangement of poles the perturbation is creating damped waves moving 
to the right for x > and to the left for x < 0. The waves propagate away from the 
origin of the perturbation at x = and are exponentially decaying with the distance 



from the perturbation. In subsection 2.1 we prove that for the holographic retarded two- 
point functions the poles do indeed fall into the first and third quadrants of the complex 
momentum plane. A typical setup with the corresponding integration contours is depicted 
in figure |2[ The imaginary part of the complex wave number can be interpreted as the 
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Figure 2: The relevant integration contours for the poles in the complexified momentum-plane. 
Figure (a) shows the contour for a; < 0, and figure (b) shows the contour for a: > 0. In order 
to obtain exponentially decaying waves travelling away from the origin of the perturbation it is 
necessary that the poles lie in the 1st and 3rd quadrants. 

inverse of an absorption length. For a given complex momentum pole g„ the right-moving 
wave has the form e~^('^^~in^) Q-ii^, The amplitude of the wave has decayed to a factor of 
1/e at a distance of A„ = l/g^. 

In the following we will be interested in computing these absorption lengths and their 
frequency dependence in the holographic dual of the = 4 supersymmetric gauge theory 
in the plasma phase. We will do this for different kinds of perturbations corresponding to 
certain gauge- invariant operators. In the gravity side we have to solve wave equations with 
purely infalling boundary conditions at the horizon just as in the calculation of quasinormal 
modes. At the boundary of AdS we have to specify the same boundary conditions that 
have been described in refs. [8,13] for the quasinormal modes. The imaginary part of these 
complex momentum wave numbers give absorption lengths characteristic of the black hole. 
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After having traveled a distance A„ a part of the wave has fallen into the black hole such 
that the amplitude is diminished by a factor of e^^. In the gauge theory dual the inverse of 
the imaginary part of the momenta can be seen as the absorption lengths for perturbations 
of the plasma by sources corresponding to gauge-invariant operators. 

Thus, we see that in both cases - relaxation times and absorption lengths - the gravity 
waves are subject to the infalling condition at the horizon. The question is simply which pa- 
rameter of the retarded Green's function is analytically continued to complex values, either 
the frequency or the momentum. To compute these complex momentum wave numbers 
one can therefore follow the same strategy that is used for the calculation of quasinormal 
frequencies, but fixing the frequency u to be real-valued instead of the momentum q. 

This switch of roles is particularly clear in the case of the AdS3/CFT2 correspondence, 
where the exact retarded Green's functions can be calculated in both sides and seen to 
match [11]. Let us consider the case of a field with conformal dimension A = 2. Then the 
retarded two-point function is 



The poles of the function determine the quasinormal frequencies ujn = ±g — i47rr(n-|- 1). 
For each quasinormal mode the dispersion relation a;„ = Un{q) is linear. Because of this 
linearity the poles can also be interpreted in a different way by writing 



where we see explicitly that the complex momentum modes lie in the first and third quad- 
rants for the right- and left-movers respectively. 

In higher dimensions the dispersion relations for the quasinormal frequencies are not 
linear and can be computed only numerically. At zero momentum, the position of large 
frequencies in the complex momentum plane has been estimated using semiclassical meth- 
ods [24], it would be interesting to extend those analysis to non-zero momentum. . Since 
the dispersion relation for the quasinormal modes is known only numerically we also have 
to resort to numerical methods to find the complex wave numbers and absorption lengths. 
The only exception is given by the hydrodynamic modes that appear for small frequency 
and wave numbers [25,26]. We will see that our numerical results are in agreement with 
the analytic dispersion relations of the hydrodynamic modes. 

2.1 Stability analysis 

We will now perform a stability analysis analogous to the one for quasinormal modes 
in [10]. We will see that the complex momentum wave numbers indeed lie in the first and 
third quadrants of the complex g-plane for positive frequencies. Note that a pole in the 
second or fourth quadrant would allow to construct outgoing waves that are exponentially 
growing with the distance from the perturbation. For the stability of the system under the 
perturbation the absence of such poles is therefore crucial. 

The time and space dependence of the field is given by simple exponentials 




(2.6) 



qn = ±[uJ + i^TTT{n+l)] , neN 



(2.7) 



^{t, x) ~ e 



—iLdt iqx 
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where q := + iq^ . We will distinguish between the cases x > and x < 0. 

For X > stability demands an exponentially decaying wave and therefore signg^ = 
+1. We further demand that the wave is outgoing from the origin of the perturbation 
which demands sign a; = signg^. Taking these two facts together amounts to the condition 

sign = sign/ . (2.8) 

Doing the same analysis for x < 0, one finds that the perturbation moves away to the 
left if signo; = — signg^, whereas the stability condition is now signf;^ = —1. This again 
amounts to eq. (|2.8|) . 

We want to prove now that in the gravity dual the complex momentum modes of the 
black hole follow indeed the rule given by equation ( |2.8D . We consider a minimally coupled 
scalar $ in AdSs x with mass m. The line element of the AdS black hole with planar 
horizon is 

d4dS = J (-/(Odt^ + dx2) + dr^ , (2.9) 

where /(r) = 1 — (rn/T)^. The temperature is given through rn = irR'^T. We will 
use in the following the coordinate z = r}i/r and rescale time and space coordinates by 
rHi?~^(t,x) 1-^ (i,x). The boundary is now located at z = and the horizon at z = 1. 
The equation of motion for a minimally coupled scalar ^{t,z,x) = exp{—iujt + iqx) ^{z) 
of mass m is 

We can further split ^{z) = a{z)y(z), in order to find an equation for y{z) that is 
Schrodinger-like in a 'tortoise' coordinate defined through 

dz* = -^ (dl+uj^ -V{z,))y{z,)=0 , (2.11) 

provided that a{z) fulfils 



a'(z) 3 



a{z) 2z ■ 

In the z coordinate the Schrodinger potential reads 



(2.12) 



V{z) = li^(^l5 + + ^ g^4^ ._ y^(^) ^ j^g(^2) ^ ■ Im(^2)y.(^) ^ (2.13) 

where we have separated it into its real and imaginary parts. In the coordinate the 
horizon lies at z* +00 and the potential vanishes there, so the wavefunction can be 
described as a superposition of plane waves. The infalling boundary condition corresponds 
to setting y(z^,) = e**^^* xi^*) with x(+oo) = const. Thus we find 

(dl + 2iu;d,^ - Vo{z.) - Re(g2) /(z,) - ilm(g2) /(z,)) ^(z,) = . (2.14) 
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If we multiply by the conjugate xi^*) pick out the imaginary part of the equation we 
obtain 

ixdlx - Xdlx) +u^O,,\x\' - lm{q^)f{z,)\xf = , (2.15) 

Now we integrate this equation between the boundary {z.^: = z^) and the horizon {z^ = +00). 
Upon a partial integration the derivative terms cancel each other: xi^*) vanishes at the 
boundary due to the Dirichlet boundary condition we impose there and at the horizon the 
derivative vanishes dz,xi+^) — 0- The remaining terms in equation ( p. 15 ) amount to 



uj \x{z = l)V= Im(g^) / dz, f{z,)\x{z.)V = M^") / dz \x{z)\' , (2.16) 



Jo 



The integral on the right hand side is positive definite, which then implies that signcj = 
signlm((7^) = sign(g^q'^) which is precisely the stability condition ( |2.g| ). 

There is a further stability condition involving the properties of the potential. When 
a; = 0, we have the condition that Im(g^) = 0, so either = ov = 0. Consider now the 
real part of the Schrodinger equation ( p. 14 ). After multiplying by x{z*) and integrating 



between the boundary and the horizon we find 

dz, (|a,,x(^*) + iuJX{z.)\^ + (Voiz,) - cv^)\x{z*)\^) = - Mq^) / dz \x{z)\^ , (2.17) 

Jo 

Clearly, if Vo{z^) > between the boundary and the horizon, then, at w = 0, Re(g^) < 
and we will have q^ = and q^ 7^ 0. On the other hand, if the potential is negative 
on some region, then there could be solutions with Re(q'^) > or equivalently q^ ^ 
and q^ = 0. Considering four-dimensional Minkowski slices of AdSs, these modes can be 
regarded as tachyonic instabilities of negative mass squared = —{q^)"^- Notice that 
with our choice the boundary conditions ~ g*'^^* and fixing lo to be real, this condition 
actually refers to the presence of 'negative energy' modes in the scattering spectrum, so 
only when the potential is negative at the horizon this kind of instabilities could appear. 
Other instabilities associated to the presence of bound states could be present, see the 
appendix ^ for a discussion. 

If Re(g^) > instabilities are present in the bulk theory, the gauge correlation functions 
associated to the dual operators will show an oscillatory behavior at large separations, as 
opposed to vanishing, indicating that the plasma is actually out of equilibrium. From the 
point of view of the effective three-dimensional theory, instabilities will appear as tachyonic 
states in the spectrum. 



3. Absorption lengths: numerical analysis 
3.1 Scalar operators 

As a first example we want to compute the absorption lengths of a scalar operator ©(t, x) 
of conformal dimension A = 4. We choose this particular case because it is the simplest 
setup we can use to illustrate the method, since the dual supergravity field corresponds to 
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a minimally coupled, massless scalar. A possible example is given by O = tvF^, that maps 
to the dilaton in the holographic dual. 

Consider the retarded two-point correlation function in the theory at temperature T 

Gn{t-t',^-y) = -ie{t-t') {[Oit,^),Oit',y)]) . (3.1) 

At large distances |x — y| ^ T~^, the Green's function decays exponentially due to thermal 
screening. As we have explained, this behavior is determined by a set of discrete lengths 
that in linear response theory describe the absorption of out-of-equilibrium perturbations. 
For the theory in equilibrium they are identified with correlation lengths in the plasma, 
as proposed in ref. [20]. The squared inverse of the zero-frequency correlation lengths can 
also be regarded as the glueball masses ^ of a three-dimensional effective theory in the high- 
temperature limit [17-19]. Via AdS/CFT correspondence we can reduce this complicated 
non-perturbative problem in the gauge theory to finding the complex momenta that allow 
the dilaton fluctuations to obey infalling boundary conditions on the horizon and Dirichlet 
ones on the boundary. In this example, and in the other cases we consider in this paper, the 
equations of motion can be reduced to Heun equations that we can solve using semi-analytic 
methods. 

The equation of motion for this field was already derived in subsection 2A . Throughout 
the paper we use dimensionless frequency and momentum. In order to recover the dimen- 
sionful quantities it is enough to make the substitution (a;,q) i-^ iTT{uj,q). Changing 
coordinates from z to x = 1 — z^, the equation now reads 

^ l+_(l-x)^ ^ / ^ ^ ^^^,(x) = 0. (3.2) 

x{l-x){2-x) V4a;2(l-x)(2-x)2 4x{l - x){2 - x) J ^ ^ 

This equation has four regular singular points at x = 0, l,2,oo, with characteristic 
exponents 

{0; -iw/4, +icj/4} , {1;0,2} , {2; -w/4, +a;/4} , {cx);0,0}. 

Therefore, we can transform it into a Heun equation and we can follow the analysis de- 
scribed in [12]. To compute the complex wave numbers we simply have to analytically 
continue the momentum instead of the frequency. It is interesting to observe that none of 
the characteristic exponents at the singular points depend on the momenta. We factorize 
$(x) as 

<l>{x) = x-'^l^ (1 - xf (2 - xyl^ y(x) , (3.3) 

which allows us to write the equation of motion in the standard form of a Heun equation 
for y{z) 



with parameters 
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a = /? = 2--(l + i) , Q = 4 + ^-(l + 7z)--(2-i)— , (3.5) 
7 = 1- i|, 5 = 3, e = l- |- (3.6) 



"in this particular example we are considering J^^ — 0^^ glueballs. 
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All the other perturbations we will consider in this paper can be transformed to Heun 
equations in a similar way. In [13] the same perturbations have been studied in order 
to solve for the quasinormal modes. The solution yo{x) corresponding to the infalling 
boundary conditions at the horizon x = is a linear combination of local solutions yi^A, 
yi^B at X = 1 

yo{x) = Ayi^A{x) + Byi^B{x) , (3.7) 
where yi^A is analytic at x = 1. The retarded Green's function turns out to be proportional 

*° A 

Gr oc - , (3.8) 

and the poles correspond therefore to the solutions with i? = 0, i.e. solutions that are 
analytic in the interval x S [0,1]. These boundary conditions determine a discrete set of 
complex momentum eigenvalues if we fix the frequency lo to real values. 

We can find local solutions using the Frobenius method close to the singularities. In 
the cases under consideration, a solution with boundary condition y(0) = const, will be 
a superposition of solutions with exponents 1 — (5 < and close to the AdS boundary 
(x = 1). The boundary condition y(l) = const, can be satisfied only for a discrete set of 
frequencies u: or momenta q. These values can be computed imposing matching conditions 
at some intermediate point for the Frobenius series although we will need a large number of 
terms and the convergence gets worse for higher frequencies. There is an alternative method 
based on the improved convergence of the solutions. Normal solutions are convergent for 
|x| < 1, but for some values of the parameters the solutions can converge for |x| < 2. 
This condition of extended convergence boils down to a transcendental equation for the 
frequencies or momenta in the form of a continued fraction (see [12, 27] for more details) 
using Pincherle's theorem on the existence of minimal solutions to three term recursion 
relations. 

The coefficients of the Frobenius series at x = should satisfy the recursion relation 
an+2+ An{uj,q)an+i+ Bn{LO,q)an = Q , n>0, (3.9) 



where 



A.(..,)=- '"^if,-^'^!'?r,»^'^ . (3.10) 

2{n + 2){n + 1 + 7) 
R ^ ^ {n + a){n + l3) 

Bn{^,q) = ^, , , , ^ , (3.11) 

2{n + 2)(n + 1 + 7) 



and oq = 1, fli = (5/27. Then, using the recursive definition 

^ _ Qn+i _ Bn{uj-,q) 

an An{uj,q) +r„+i 



(3.12) 



Pincherle's theorem states that a minimal solution to the three term recursion relation 
( 3.10 ) exists if and only if the continued fraction on the right hand side in ( |3.12| ) converges. 



Moreover, in this case it converges precisely to a„+i/a„. In [12,13] it was pointed out 
that the minimal solution corresponds precisely to a solution of the Heun equation that is 
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analytic at x = 1 therefore fulfilling the correct boundary conditions. Choosing n = we 
find 

(3.13) 



Q 

27 



and computing recursively gives a transcendental equation for q. Using this formula, we 
can compute numerically the complex momentum modes with high precision. In order to 
do that, we cut the fraction at a large value n = = 100 and use the asymptotic value 
r„ = 1/2 — (2 + a;)/4n*. It is important to realize that Pincherle's theorem applies only if 
we are dealing with genuine three term recursion relations. Sometimes it can happen that 
the recursion relation involves three terms only from a certain value of n = ni on. This 
happens for example if either a = or /3 = when Bq = 0. In such a case one has to use 
(3.12) with n = ni. We will see that we are faced with this in the cases of the longitudinal 
vector field perturbations and of the shear mode perturbations at a; = 0. Since q = in 
both cases it is sufficient to take ni = 1 and use 

Q2 + 3Qj - 2a/37 + 2Q5 + Qe 



ri 



4Q + AQ-f 



(3.14) 



instead of ( p. 13 ). 



qj versus qu for the first five complex momentum eigenvalues 
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Figure 3: We have traced the locations of the five lowest momentum eigenvalues in the complex 
g-plane for different frequencies as a function of the frequency out to w = 50. The momentum 
eigenvalues vary continuously with the frequency and lie on the analogues of Regge trajectories. 

We have numerically computed the complex momentum eigenvalues using this method. 
The results for the scalar field perturbations are shown in figure The real and imaginary 
parts of the five lowest complex momentum are plotted as a function of the frequency. 
The real parts start out at = for zero frequency. The imaginary parts start out 
at a finite value at u = 0, develop a shoulder that is more pronounced for the higher 
modes and then fall off rather fast until they enter a regime of slow decrease for large 
frequencies. Numerically we found that the lowest mode becomes almost constant at 
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large frequencies with ~ 0.83 at w = 100. Also the higher modes flatten out for high 
frequencies. As expected, higher frequencies can penetrate farther into the plasma. It is an 
interesting question if the plasma becomes transparent for some high but finite frequency, 
if transparency is reached only in the limit u; = oo or if the absorption length stays finite. 
Unfortunately our algorithm does not allow us to explore this asymptotic regime. We can 
speculate however using the underlying conformal invariance of the AA = 4 theory. Since for 
high frequencies the temperature is less and less important we expect that the absorption 
length diverges as a; — > oo, i.e. q\u} = oo) = 0. A finite absorption length would point 
to an underlying scale in the theory. On the other hand, if the plasma were to become 
transparent at some finite value of w, we would expect that to happen at a scale that is 
set by the temperature. However, our numerical results show finite absorption lengths for 
much higher frequencies. 



3.1.1 Glueball masses 



Of particular interest are the absorption lengths in the static limit 
cj ^ 0. In this case we will refer to the absorption length as the 
screening length. The equation ( ^^ ) with u; = has been studied 
before in [17, 18]. There the interpretation of the eigenvalues in 
the momentum with < Q was as masses of glueballs in the three 
dimensional theory that is obtained by reduction on the thermal 
circle in the Euclidean section of the AdS black hole. The glueball 
masses can be calculated as the discrete eigenvalues = —q^. Our 
numerical results at = for the first ten modes are compiled in 
table H and are in good agreement with results given by refs. [17, 18]. 

It is important to see if the eigenfunctions correspond to the 
wave functions of the glueballs too. In [19] the authors observe that 
for all the glueball masses the correct boundary conditions corre- 
spond to demanding analyticity of the wave function at the horizon 
and the boundary. These are precisely the same boundary condi- 
tions that emerge in our case at w = 0. Therefore the screening lengths for static fields 
corresponds precisely to the glueball masses computed earlier in [17-19]. 



n 




1 


-11.5877 


2 


-34.5270 


3 


-68.9750 


4 


-114.9104 


5 


-172.3312 


6 


-241.2366 


7 


-321.6265 


8 


-413.5009 


9 


-516.8597 


10 


-631.7028 



Table 1: The first ten 
glueballs of the scalar 
mode. 



3.2 Global currents 

In the AA = 4 theory, the global currents associated to R-charges map to mixed components 
of the and AdSs metrics, that can be seen as graviphotons after dimensional reduction 
to AdSs. In general, any global symmetry in the field theory will map to a local gauge 
symmetry in the holographic dual. Then, to find the poles of the retarded Green functions 
in the plasma 

G^,{t-t',i^-y) = -ie{t-t') {[J^{t,iL),Mt',y)]) , (3.15) 
we have to compute the complex momentum eigenvalues for vector fields in the AdS black 
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hole.^ We will see that there are two decoupled sectors, corresponding to transverse and 
longitudinal channels. The reason is that temperature breaks four-dimensional Lorentz 
symmetry to three-dimensional rotational symmetry. In the glueball language, the zero- 
frequency masses correspond to J^*^ = 1 and ^ states. However the states arising 
from the vector fields are also charged under the global R-symmetry and therefore do not 
form part of the superselection sector that constitutes QCD3. For simplicity we will refer to 
these states also as glueballs. The longitudinal channel is special because it also describes 
the diffusion of the conserved charge through the plasma, that does not appear as a glueball 
state in the three-dimensional theory because the residue of the diffusion mode vanishes in 
the zero- frequency limit. We will show that the diffusion pole is also captured by complex 
momentum eigenvalues. 

We can compute the complex momentum eigenvalues corresponding to a vector field 
in the AdS-Schwarzschild background in an analogous way to the scalar field case. The 
equations of motion for such a field are given by the Maxwell equations 



=5. [y^^'^Vi^p.] = 0, (3.16) 



where Ff^^, = d^Ay — dyA^. We can choose the = gauge in the metric (|2.9|) and expand 
in plane waves A^. Separating the vector field in longitudinal and transverse components, 
Al{z) = ^ ■ A{z) and q • Ar(z) = 0, the equations of motion are 



Q = uA'^{z) + f{z)qA'L{z), (3.17) 



1 . 1 / ... 2 



= A'^{z) - -A'^(z) - —-iuqALiz) + q'A^{z) , (3.18) 
z f{z) V J 

= AUz) + Al(z) + jl^(u;qAo{z) + u;'AUz)) , (3.19) 











z 




f{z) 





= AJ(.) + 7^ Ak.) + [J^^ - Ti;) j At(^) . (3.20) 

The first three equations are not independent so we can use the first one to write decoupled 
equations for Aq{z) and A'j^{z). Notice that there is a choice of gauge invariant variables 
El = qAo + ivAl and E-p = loAt that describe the diffusive and transverse channel 
respectively [8]. However, the spectrum of complex momentum values (equivalently of 
quasinormal modes) is gauge invariant, so it should not matter if we choose to work with 
gauge components, that obey simpler equations. Since the invariant quantity is El, this 
means that Aq and Al should have the same spectrum, as the constraint ( 3.17| ) points out. 



Then, the relevant equations for the temporal, longitudinal and transverse components 



^We are assuming that the total charge in the equihbrium state vanishes, so there are no chemical 
potentials. 
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of the vector field read 

= A^^ + f 4 - - Vt + - ^1 At(z) . (3.23) 



wliere we define Vo{z) = Aq(z), Vl{z) = A'^(z). In the x = 1 — z^ coordinate and for 
a suitable factorization of each component, the equations above can be written as Heun 
equations: 

Temporal. The critical exponents at the singularities are 

1 3 



Vq{x) = x-'^^/^ix - 1)^/2(2. _ 2)-^/4y(x) . (3.24) 



Longitudinal. The critical exponents at the singularities are 

(0;-l-^^,-l + ^9 , (l;^,^), (2;-l-^,-l + 9 , (oo;^-A^ + V5 

Vl{x) = x-'^~'^/\x - l)^l^{x - 2)-^-^'S{x) . (3.25) 
In both cases we find the same parameters for the Heun equation 



a = --(l + i), /3 = 2-^-(l + i)j , Q = :^-(l + 3i)--(2-i) ^ 
7 = 1- 5 = 1, e = l-|- (3.26) 

Notice that the boundary conditions for Vl{x) are not infalling ones. They are determined 
by the constraint ^.17 . 

Transverse. The critical exponents at the singularities are 

0;-^^,.^) , (i;0,i) , (2;-^,^) , (oo;0,i) . 

At{x) = x-'^l\x - \){x - 2)-^l^y{x) . (3.27) 



a = l--(l + z), /3 = 2- (^-(l + i)j , g = :^ + 2-(l + 5i)--(2-i) ^ 

7 = 1- i|, <5 = 2, e = l-|- (3.28) 

As we had anticipated, the temporal and longitudinal components have the same 
spectrum, since they obey the same Heun equation, although this was not evident in 
equations ( |3.21[ ) and ( |3.22| ). 
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The results are shown in figures |5| and |^. The real and imaginary parts of the five 
lowest complex momentum are plotted as a function of the frequency. The behavior is 
similar to the one found for the scalar operator. The imaginary parts start out at a finite 
value at a; = 0, develop a shoulder that is more pronounced for the higher modes and then 
fall off rather fast until they enter a regime of slow decrease for large frequencies. The real 
parts start out the = for zero frequency. 

So far, we have described the absorption of R-current excitations in the plasma. How- 
ever, a conserved global charge cannot be dissipated, it is spread out by the slow process of 
diffusion. This is described in the hydrodynamic regime w, g <C T by a diffusion pole [25] 
(units restored) 

- = -^. (3.29) 

In our analysis of complex wave numbers we are able to see numerically this mode {q = 
(1 + i)y/uj with our conventions) that fits nicely with the analytic prediction in the hydro- 
dynamical regime Fig. |8|. 

3.2.1 Glueball masses 

In the zero frequency limit, the absorption lengths can be interpreted as the inverse glueball 
masses of an effective three-dimensional theory. Note however that these states do not 
lie in the superselection sector that constitutes the holographic dual of QCD3! For the 
longitudinal channel we have to take into account that a{uj = 0) = so we have to use the 
modified recursion relation starting at n = 1 ( 3.14| ). It turns out that the glueball masses 



of the longitudinal channel coincide with the ones found for the scalar operator, table |T|. 
Indeed, at u; = we can transform the Heun equation with parameters (|3.5| ) into the Heun 
equation with parameters ( 3.26 ). First, make the coordinate transformation x — > fEj j that 
shuffles the singular points 2 0, 1 00. Then, the redefinition y{x) {x — l)'^y{x) 
(explained in the appendix) shows that both equations are equivalent. Notice that the 
solutions that are analytic in [0, 1] in the transformed equation correspond to solutions that 
are analytic in [2, 00] in the original equation, and not to the physical modes. However, 
such solutions can be generated from the physical ones by conformal transformations on 
the two sphere^, so both types appear for the same values of the parameters. Notice that 
the both solutions have a similar analytic structure, the only singularity is a branch cut 
joining two of the singular points. Also the fact that the auxiliary parameters Q of both 
equations are the same for the particular cases we are considering, allows an immediate 
identification of the complex momentum numbers. 

The transverse channel has different spectrum, whose first modes are in table Al- 
though the glueballs associated to vector fields have non-zero R-charge, and are usually 
not considered, our computation shows that the lightest three-dimensional state and hence, 
the longest correlation length, belongs to this class. ^ 

3.3 Stress-energy tensor 



See [30] for an exhaustive list of Heun solutions and their relations. 
^This state is even lighter than the lightest QCD3 glueball listed in [19]. 
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The stress-energy tensor of the gauge theory encodes important 
dynamical and thermodynamical properties of the plasma. Correla- 
tion functions of the stress-energy tensor 

G^,,p^{t-t',x-y) = -ie{t-t') {[T^^{t,^),Tp„{t',y)]) , (3.30) 

are related to perturbations of the metric that leave the factor 
invariant. Therefore, we want to introduce a small fluctuation of 
the four-dimensional part of the metric g^i, = g^^ + h^y . 

In the gauge theory, the breaking of Lorentz symmetry to ro- 
tational symmetry by temperature splits the Green's functions in 
transverse, shear and sound channels, that in the zero frequency 
limit contain the J^'~" = 2++, and 0^^ glueball spectrum. This 
is reflected in the gravity dual, where the perturbations fall into 
three different classes with decoupled field equations [8,25,26]. The 
associated spin to each of these channels is also 2, 1 and 0, so we 
will refer to them also as tensor, vector and scalar. 

In the shear and sound channels there are also hydrodynamical modes that describe 
the diffusion of conserved momentum and the propagation of sound. We will not study the 
sound channel, but we will show that complex momentum modes also capture the shear 
pole. 

We will work with gauge-invariant variables, following [23]. There, the authors consider 
general metrics of the form 



n 




1 


-5.1313 


2 


-22.4816 


3 


-51.2098 


4 


-91.4106 


5 


-143.0926 


6 


-206.2577 


7 


-280.9066 


8 


-367.0395 


9 


-464.6566 


10 


-573.7580 



Table 2: The first ten 
glueballs of the trans- 
verse mode. 



ds^ = -F(r) dt^ + -— + r'^da. 

F(r) 



(3.31) 



where da^ corresponds to a metric of a n-dimensional space of constant sectional curvature 
K = 0,±1, and 

Hr) = K-^-Xr'. (3.32) 



^n—l 

In our case, K = 0, n = 3, X = —1 and M = 

The Einstein equations are decomposed in tensor, vector and scalar components rel- 
ative to the three-dimensional metric. It is thus possible to define three different gauge- 
invariant quantities to which we can associate a Schrodinger-like equation of motion [23]. 
In the z coordinate it reads 



dz 



+ Vi{z) 



I = {T,V,S} , 



(3.33) 



where for each perturbation we will have a different potential. Rewriting the Schrodinger 
equation by shifting {LO,q) ru{i^,q), the potentials V are given by 



Vriz) 
Vv{z) = 
Vs{z) = 



4z2 

m 

4^2 

M. 

4^2 



(15 + 4^2^2^9^4) ^ 



(3 + 4^2^2 



27z' 



(3.34) 



(1 + 6g-2z2)2 



2 4 8 ' 

1 + 4:q^z'^ + Qz"^ + 156z^ - 108^ + 540^ + 324^ 

Q Q Q , 
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for tensor, vector and scalar perturbations, respectively. 

By making the change of variable x = 1 — z^, the equations for tensor and vector 
perturbations lead to a Heun equation. For scalar perturbations the situation is not so 
simple, and it requires a separate analysis that we leave for future work, so in the following 
we will be concerned only with tensor and vector perturbations. 

Tensor perturbations. The characteristic exponents are 

uj ,L0\ ( ^ 3 5\ / uj uj\ ( 3 3 



^^(x) = x~'^l\x - \fl^{x - 2ylS(x) . (3.35) 

a/3= (-(l + i)-2j , Q = ^ + 4-(2-i)(^(-l + 3i)- + — 
7 = l- z|, 5 = 3, e=l-|- (3.36) 

Vector perturbations. The characteristic exponents are 

/ u) ( 13\ ( ijj LO\ ( 39 

^y{x) = x-'^l^{x - \fl^(x - 2)-^/^2/(x) . (3.37) 

a/? = ^(l + 0(^(l + i)-3) , Q = ^-(l + 5i)^-(2-z)^^ , 

7 = l-i|, 5 = 2, e = l-|- (3.38) 

Notice that tensor fluctuations obey the same equations as a massless scalar field, so the 
first modes of the spectrum are plotted in figure ^. The Heun equation we have for the 
vector perturbations goes over to the one the authors in [13] found for the shear mode after 
the transformation described in the appendix. 

As we have commented above, vector fluctuations correspond to the shear channel of 
the gauge theory. This channel is associated to the momentum of the plasma, that as a 
conserved quantity is not absorbed but diffused. In the hydrodynamical limit it is possible 
to find an analytic expression for the diffusion pole [25] 

uj = -i . (3.39) 



We find good numerical agreement for this mode g = (1 + i)v2a;, as can be seen in fig. ^. 

The results for the shear mode are shown in figure |^ The real and imaginary parts 
of the five lowest complex momentum are plotted as a function of the frequency. Again, 
we find a similar behavior to scalar and vector modes. The imaginary parts start out at 
a finite value at a; = 0, develop a shoulder that is more pronounced for the higher modes 
and then fall off rather fast until they enter a regime of slow decrease for large frequencies. 
The real parts start out the (7^ = for zero frequency. 
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3.3.1 Glueball masses 



We can find the glueball spectrum of the effective three-dimensional theory by taking the 
static limit a; = 0. Again we have to use the recursion relation starting n = 1 ( ^.14 ) since 
a{ijj = 0) = 0. The results for the shear channel are compiled in table ^. The glueball 
spectrum for neutral glueballs has been computed using a similar supergravity approach 
in [19]. The numbers we find differ actually somewhat from the ones quoted in [19] for the 
glueballs. We attribute this to the different numerical methods that have been used 
to obtain them. 



4. Conclusions and Outlook 



n 


9n 


1 


-18.6758 


2 


-47.4951 


3 


-87.7228 


4 


-139.4167 


5 


-202.5882 


6 


-277.2408 


7 


-363.3762 


8 


-460.9949 


9 


-570.0974 


10 


-690.6838 



Table 3: The first ten 
glueballs of the shear 
mode. 



We have established a relation between solutions to linearized field 
equations with complex momenta in an AdS-black hole background 
and the absorption lengths of a conformal gauge theory in a plasma 
phase. We have explicitly studied some simple examples correspond- 
ing to scalar, vector and metric fluctuations. Due to conformal 
symmetry, all absorption lengths scale simply with T~^. At zero 
frequency we find agreement with previous computations of the ef- 
fective three-dimensional glueball spectrum [17-19]. However, we 
prefer in this paper to interpret our results as screening lengths for 
static fields. This interpretation has also recently and independently 
been proposed in ref. [20]. 

Furthermore, we have computed the dependence of the absorp- 
tion length on the frequency. The results for the first modes are 
compiled in figures ^, ^, ^ and 0. In all the cases, the plasma is 

less absorptive for higher frequencies. The complex wave numbers also capture the hydro- 
dynamical behavior for R-charge and momentum diffusion. Our numerical results are in 
agreement with the simple analytic continuation of the dispersion relation for the hydro- 
dynamic modes. This is shown in figures ^ and 

One of the interesting results of our study is that the longest screening length (the 
lightest "glueball" mass in the dimensionally reduced theory) corresponds to a state with 
non- vanishing R-charge. Such a state does not belong to the spectrum of the QCD3 theory, 
i.e. the mass gap of the effective three dimensional theory is not the one of QCD3! Glueball 
masses play an important role in the determination of the Debye screening length. Here 
one studies the glueball exchange between open strings in the AdS black hole background. 
As has been pointed out in [20] the mass gap by itself is not important for the Debye 
screening, because only specific operators can couple to the open string. Since these open 
strings are R-charge neutral, the low mass states with non zero R-charge do not couple to 
the string. However, the string configuration one considers usually has its endpoints fixed 
on one point on the S^ and it is also possible to consider strings that end on different points 
on the S^. In such a situation the light non-zero R-charge states might become relevant 
and could modify the result for the screening length. 
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In this paper we have only studied the cases that can be reduced to Heun equations 
and allow the application of the efficient continued fraction approach to the calculation 
of the complex momentum eigenvalues. It would certainly be interesting to extend the 
present investigations to the cases that cannot be reduced to Heun equations. In these 
cases one has to resort to elementary method of Frobenius expansions and this slows down 
the numerical calculation considerably. Nevertheless we think that this is an interesting 
problem especially in view of the comparison to the glueball mass calculations. 

Another rather intersting point is the question wether the absorption length diverges 
in the limit of infinite frequency or wether it stays finite. Unfortunately so far we know 
only about numerical methods to evaluate the absorption lengths. 

A related problem is the calculation of the the absorption lengths in non-conformal 
holographic theories. Due to the presence of an underlying scale the dependence on the 
frequency is likely to show a more complicated pattern than the one we have found for 
the conformal case in this paper. It will also be of high interest to compute absorption 
lengths for the meson states that appear in theories with D7-brane embeddings in the 
AdS black hole using the same methods that have been employed in the study of meson 
quasinormal modes in [28]. In [29] it has recently been emphasized that instabilities arise 
for near critical black hole embeddings. Such instabilities show up as quasinormal modes 
with positive imaginary part. As we have seen, similar instabilities can also arise in the 
study of the absorption lengths. Since the instabilities in the screening lengths arise at 
a; = and for real values of it might be much easier to search for these instead of 
unstable quasinormal modes! 

We hope to make progress on these questions in future research. 
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A. Effective potentials 

In section 2.1 we have presented the stabihty analysis for a scalar field but it can be gener- 
alized for any field component (p{z) satisfying a decoupled linear second order differential 
equation 

^"{z) + Ai{z) <f'{z) + Ao{z) ip{z) + B{zfuj^^{z) = . (A.l) 
Factorizing (p{z) = cr{z) (j){z) and normalizing the (p" term 

0" + (^2^ + Ai^ (J)' + (^Ao + Ai^ + ^'^cl) + B{zfuj^ / = . (A.2) 

We now change B{z)dz = dz* and divide by B{z)'^ 

This expression becomes a Schrodinger equation when a satisfies 

Then, the same stabihty arguments can be applied with the proper identification of the 
potential 

We will now apply this to the other equations under consideration in this paper. 

• transverse vector components 

V, = M{3 + 5z^ + 4q^z-') (A.6) 

• longitudinal and temporal vector components 

Vo = -^{l + 7z'-4qh') (A.7) 

• gravitational vector perturbation (shear mode) 

l^o = ^(3-27/ + 4gV) (A.8) 

Due to the underlying analyticity of the solution of the corresponding Heun equation 
ah fields, At,Vl,Vo,^v, fulfill the boundary conditions leading to ( p^ ) ( |2T^ ). The 
effective potential is positive in the case of the transverse Vector fields. For the longitudinal 
and temporal vector field components it is negative and therefore the stability argument 
presented in section 2A does not apply. We note that the asymptotic behavior at the 
boundary is the same as that of a scalar field saturating the Breitenlohner-Freedman bound. 
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We take this as an indication for stability, in the original analysis in AdS a positive energy 
condition is satisfied even for fields with negative potential [31]. The asymptotic behavior 
of the fields is restricted by the condition of having a well-defined conserved energy. In turn, 
the positive contribution of the kinetic energy always overcomes the negative contribution 
from the potential. A formal analysis [32] can be applied that shows the stability of vector 
perturbations. The 'Hamiltonian' operator H = —d1^ + Vq must be positive definite over 
the set of normalizable solutions 



dz,x*Hx > . (A.9) 



We can rewrite (A.9) as 

/•oo 

- [x*Dpx]T^ + / dz, {\Dpx\^ + Vp\x\^) > , (A.IO) 

where we have introduced an auxiliary function p, such that Dp = dz, + p and 

Vp = Vo + dz,p- . (A.ll) 

A convenient election that makes > for vector fiuctuations \s p = —f{z)/2z. We 
can easily see that there is no contribution from the boundary term at the horizon, since 
x(oo) const, and p{co) 0. Therefore, we are left with the condition 

Close to the boundary, Vq ~ — l/4z^, so the solution is a combination of Bessel functions 
X ~ a^/zJo{u)z) + b^/zYQ{uJz). Then, the condition (A. 12) satisfied when 6 = 0, that is 



equivalent to choose the normalizable solution at the boundary. 

The effective potential of the shear mode is also interesting. It is negative close to 
the horizon. Equation ( 2.17| ) shows that this is a necessary requirement for existence of 



the hydrodynamic shear mode with Re(g^) = 0. We would expect that if the potential 
is deep enough, instabilities will appear. This is in agreement with other analysis that 
exhibit a negative well in the interior. Purely imaginary frequencies have been found in 
the study of electromagnetic and gravitational perturbations in global AdS [33,34]. In the 
extremal limit, the frequencies seem to reach the real axis at a; = 0, and the geometry was 
conjectured to be marginally unstable. Recent works also suggest that instabilities of D7 
probe branes in AdS appear when a quasinormal mode cross the real axis at a; = [28,29]. 

B. Changing parameters in a Heun equation 

In this appendix we show how to map a given Heun equation with given parameters into 
another Heun equation for a different function with a different set of parameters. This will 
allow in some cases to avoid the problem with "fake" modes. 
Let us start with a Heun equation for y{x) 

" + ("^ + ^ + ^ '^^^ ~ ^ ( ) ~ (B 1) 

\x x — 1 x — 2) x{x — l){x — 2) ' 
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where the parameters are subject to the condition a + P+1 = j + S + e. The characteristic 
exponents at the AdS boundary (x = 1) are in general 

{l;0,l-<5} , 

for y{x). In the cases where 5 G {0 U Z~}, the first solution is logarithmic, but the 
logarithm might accidentally vanish at the "false frequencies" 

y{x) = B[l + {l-x) + ... + h{l- xf-^^ log(l -x) + ...]+ A[{1 - xf-^^ + ...], (B.2) 

where h = 0. According to [13], the Green's function is proportional to the ratio A/B, and 
the false frequencies arc the ones for which accidentally h = 0. 

Now we would like to change the function y(x) such that we find a related Heun 
equation with different parameters. Let us define 

y(x) = (l-x)^y(x) . (B.3) 

This allows us to find a Heun equation for Y{x), provided g = 1 — 6, and where the new 
set is 

a = a + g = a + l — 6, 

Q = Q + 2jQ = Q + 27(1 - 5) , (B.4) 
6 = 6 + 2q = 2-6, 
1 = 1, e = e . 

The interesting thing about this shift is that we can find a positive 6 G Z+ when in the 
original Heun equation we encounter fake frequencies. This always eliminates the false 
frequencies since the second solution is never analytic but goes like (1 — x)^^~^\ which is 
a negative (or zero) exponent for 6 G Z+. Now, the Probenius solution is 

Y{x) =B[{1- x)-^ + ... + h{l- x)^^-^-^^ log(l -x) + ..] +A[{1 - x)^^-^-^^ + ...]. (B.5) 

The recursion algorithm of Leaver [27] and its adaption to the Heun equation by Starinets 
[12] computes when the solution of the Heun equation is analytic at x = 1. Now the solution 
that goes with the coefficient B is never analytic, and therefore the false frequencies do not 
appear. 
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Req vs. u for the lowest five Scalar modes 
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Figure 4: Real and imaginary parts of the lowest five complex momentum eigenvalues versus the 
frequency. In the lower-right corner of the first figure we have zoomed in to show the separation 
between the five modes. 
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Req vs. lj for the lowest five Longitudinal/Temporal modes 
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Figure 5: Longitudinal perturbations of the vector field on AdS. At a; = the values coincide with 
the ones of the scalar field perturbations. For w > the shape is however different. 
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Figure 6: The complex momentum eigenvalues of the transverse vector field components. Note 
that the lowest mode gives the longest absorbtion length. The plasma is most transparent to 
transverse vector perturbations. 
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Re q vs. lo for the lowest five Shear modes 
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Figure 7: The lowest five complex momentum eigenvalues above the diffuse mode for the shear 
channel perturbations. 
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Re 5 vs. u) Imgr vs. u) 
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Figure 8: Diffusion mode. The solid line represents the numerical solution while the dotted line 
is the analytic formula from hydro dynamical analysis. 



Re 5 vs. u) Imgr vs. w 
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Figure 9: Shear mode. The solid line represents the numerical solution while the dotted line is 
the analytic formula from hydrodynamical analysis. 
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